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We validate that ofT-resonant electron transport across ultra-shoTt oligomer molecular junctions is characterised 
by a conductance which decays exponentially with length, and we discuss a method to determine the damping 
factor via the energy spectrum of a periodic structure as a function of complex wavevector. An exact mapping to 
the complex wavevector is demonstrated by first-principle-based calculations of: a) the conductance of molecular 
junctions of phenyl-ethynylene wires covalently bonded to graphitic ribbons as a function of the bridge length, 
and b) the complex-band structure of poly-phenyl-ethynylene. 



1. Introduction 

Since the pioneering theoretical proposal of 
a molecular rectifier by Aviram and Ratner 
molecular electronics has received con- 
siderable attention 0. A first key experi- 
ment [31 investigated the predictions of tun- 
nel theory. It established the exponential de- 
crease of the conductance in Langmuir-Blodgett 
[CH3(CH2)Ar_2COO]2Cd films between Al elec- 
trodes versus their thickness for chain lengths 
A''=18-22. Most recently, in an effort to illu- 
minate the mechanism of electron transport in 
relatively short mono-molecular junctions, such 
length dependence has been reported in experi- 
ments with scanning probe microscopy techniques 
|4|5l(il7j , where a small number of molecules con- 
tributes to the conductance. Structure depen- 
dent factors have been studied by looking at a 
variety of oligomers, and by systematically vary- 
ing the contact to the electrodes via the anchor 
groups or the tip load. In some cases, measured 
molecules had chain lengths ranging from one up 
to four monomer units |4I5I6| . All results point to 
through-bond off- resonant tunnelling transport. 
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In metal- molecule- metal junctions (MMMs), 
the molecular electronic states hybridise with the 
metal wavefunctions giving a finite width in the 
energy space. If no large charge transfer oc- 
curs from/to the molecule, the Fermi energy Ep 
lies within the (broadened and shifted) highest- 
occupied- and lowest-unoccupied- molecular- 
orbital (HOMO-LUMO) energies. Such an align- 
ment poses a potential barrier to traversing elec- 
trons accounting for tunnelling. In this regime, 
there is an analogy of MMMs to metal-insulator 
(or metal-semiconductor) junctions, and of the 
molecular-orbital energy broadening to the metal- 
induced gap states (MIGS) |8l()ll()lll| . 

Such states are induced by matching to the 
metal side, appear at energies within the band- 
gap of the insulator (semiconductor), and decay 
exponentially away from the junction in the non- 
metal IHj. MIGS are responsible for the overall 
finite density of states (DOS) within the HOMO- 
LUMO gap Eg for short molecular wires because 
of contributions coming from both electrodes [H]. 
For longer wires, one expects the MIGS to be 
mainly located near the interfaces on either side of 
the molecule. By looking at the positional depen- 
dence of the local DOS, such behaviour has been 
observed for Si nanowires between Al electrodes 
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|1(J| and octaneditiols between Au (111) surfaces 
|llj . In the latter case, a one-to-one correspon- 
dence between the local DOS decay parameter 
and the complex wavevector of wavefunctions in 
the forbidden energy domain of the corresponding 
one-dimensional crystal, was found. 

Such a property has been conjectured to also 
hold for the damping factor (3 of the conductance 
of oligomer molecular junctions |11I12I13] . deep in 
the tunnelling regime, where (3N ^ 1. The phys- 
ical meaning of f3 arises naturally, and the dis- 
persion relation with complex wavevector may be 
used to define an effective mass of tunnelling elec- 
trons p,l 13 . However, since most studies assume 
that (3N ^ 1, the range of validity of such tun- 
nelling characteristics is unclear. In what follows, 
we elaborate on this and demonstrate that the 
conductance decays almost always with length in 
an exponential fashion. We focus on the damp- 
ing factor of electron transport across oligomers 
which is the relevant quantifying quantity, read- 
ily deducible experimentally. In particular, we 
establish its connection to the band structure of 
the corresponding polymer when extended in the 
complex wavevector plane. Most importantly, 
this mapping holds even for very short bridges 
consisting of a handful of monomer units, as in 
experiments, and for surprisingly small values of 
the damping factor, providing an indispensable 
tool for the understanding and prediction of the 
tunnelling phenomena involved. 

2. Electron transport framework 

The low-bias and low-temperature electronic 
response of the MMMs is determined by the elas- 
tic scattering of independent electrons with en- 
ergy E for molecular wires in which electron resi- 
dent times are small compared to those of molec- 
ular vibrations 2 . Electrons see a frozen-nuclei 
configuration and, in a first approximation, the 
main contribution derives from equilibrium con- 
ditions. Polaronic effects may also be disregarded 
for molecular bridges of the size we consider JJ] , 
whereas, the response is voltage independent |15j . 
The Landauer conductance g is 

g{E) = ^ T{E), (1) 



where the transmission function T denotes the to- 
tal transition probability from incoming to outgo- 
ing states at mutually exclusive electrodes. The 
temperature constrain is relaxed by comparing 
the molecular energy scales ('~ eV) to room tem- 
perature 25meV). 
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Figure 1. Structure of diphenyl-ethynylene 
bonded to two graphitic electrodes. Open cir- 
cles denote hydrogen atoms. See Sec0]for more 
details. 

For energies far-off the molecular bridge 
resonant-spectrum and relatively long wires, the 
transmission function has been studied in both 
model systems and realistic molecular junctions 
[12116117118119) . It shows the characteristic expo- 
nential length dependence of a tunnelling process 

T(£;) = ro(^) e-'^^^)^^-!), (2) 

where N represents the number of monomers that 
constitute the molecular wire. In Sec^l we show 
that Eq|21 is of more general validity. 

The pre-exponential factor relates to the 
molecule/electrode interface; it defines the con- 
tact conductance go in units of {2e^/h). In the 
long oligomer limit, the damping factor /3 is a 
property of the bridge, related to its electronic 
structure. Based on WKB approximation a sim- 
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pie expression reads 



P{E)=2^r-^E,„ 



(3) 



where ^Jg, equals Eg^ or \E - Shomo(lumo)I- 
However, this overestimates (3 [TTtlS, , and calls 
for an elaborated approach as discussed later. 

To investigate the conductance properties of 
molecular junctions as that of FigQ] we calcu- 
late the transmission function with Green func- 
tion methods of quantum transport |2(Jj via 



T{E) = Tr [Tl{E)G-{E)Tr{E)G''(E)] . 



(4) 



Qr(a) ^ [ES -U- is the retarded (ad- 

vanced) molecular Green function dressed by a 
self-energy interaction S — + '^r- This fol- 
lows the Lowdin projection technique and takes 
into account the coupling to the 'left' (L) and 
'right' (R) electrodes. Twice the imaginary part 
of "^LiR) defines the spectral width r^jfl). Since 
we are interested in the low-voltage regime we 
use the equilibrium Green function. We have as- 
sumed a non-orthogonal basis which introduces 
the overlap matrix S, and an algebraic procedure 
with finite dimensional matrices. 

To treat the electronic structure of the oligomer 
and the electrodes, we employ a linear combina- 
tion of atomic orbitals approach parameterised by 
Density Functional Theory (DFT-TB) 21 in the 
local density approximation (LDA). In DFT-TB, 
the single-particle electronic Kohn-Sham eigen- 
states "02 of the system are expanded in a non- 
orthogonal basis set (yC^t for valence electrons lo- 
cated at the ionic positions R^, namely 



(5) 



The basis set is computed in terms of Slater- 
type orbitals in a modified atomic potential that 
slightly compresses the electron density to take 
into account that atoms are brought together in 
a condensed state. With the Ansatz of Eq[3 the 
Kohn-Sham equations for ijji are transformed to 
a set of algebraic equations 



^ ] {Hi_ii, — S^ijEi 



0. 



(6) 



The tight-binding (TB) approximation assumes 
superimposing atomic densities when calculating 
potential contributions 22 . include only 

two-centre, distance-dependent, Hamiltonian ma- 
trix elements. EqlSjis of the extended Hiickel type 
but all elements are calculated numerically with- 
out empirical parameters. Moreover, the a priori 
parameterisation of H and S provides an efhcient 
scheme for calculating the Green functions. 

3. Complex-band structure 

Assuming that interfacial effects are not domi- 
nant, any evanescent states in an underlying peri- 
odic system are understood as bulk Bloch states 
with a k vector with an imaginary component, 
i.e., k = q — iK. These states are quantified by the 
complex-band structure, the extension of the elec- 
tronic band structure to include complex Bloch 
vectors. In band structure calculations, real k in 
the first Brillouin zone are given as input to solve 
the energy eigenvalue problem. For the complex- 
band structure, the inverse problem is defined. 
All k vectors associated with a real energy E are 
sought, which are in general complex. The ana- 
lytic properties of the energy spectrum of crystals 
as function of a complex variable have been ini- 
tially studied by Kohn 23 . 

In any case, EqEl in a periodic potential is 
solved, which within our TB approximation reads 



ika I TT ^ — ika\iE 



„-ika\iE 



(7) 



The matrices H^j and Sij contain the Hamilto- 
nian and overlap elements, respectively, between 
the zth and jth monomer units; a is the lattice 
constant. Furthermore, we employ a useful algo- 
rithm based on transforming Eq[7|to a generalised 
eigenvalue problem for A = e*'°", yielding 



I 

X ( "{H.m.m+l ^ ESm.m+l) 



(8) 



The procedure resembles finding the eigenvalues 
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of the transfer matrix in a non-orthogonal lo- 
cahsed basis set but it overcomes problems arising 
from possible singular matrices. This formal sim- 
ilarity also justifies a relation of f3 to the complex 
wavevector. We solve numerically Eq|Sl for any 
real E and classify its eigenvalues. For |A| — 1, k 
is real; otherwise k is complex and describes ex- 
ponentially decaying or growing solutions, which 
are of particular interest in the tunnelling regime. 

4. Prototype system 

The molecular junctions under investigation 
consist of oligo-phenyl-ethynylene (OPE) wires of 
small varying length connected to two graphene 
electrodes in a ribbon shape. Before calculat- 
ing the conductance of systems as that depicted 
in Fig^ all the formed molecular junctions are 
optimised. More precisely, the non-shaded area 
is dispatched and the consequent edges are re- 
placed by single hydrogens. To simulate the 
bulk electrodes, a constrained atomic relaxation 
is performed. Graphene atoms away from the 
OPE/ribbon interfacial bond by approximately 
3 X hexagonal lattice constant remain immovable. 
The shaded-area represents infinite graphitic rib- 
bons whose structure is assumed not to differ 
from that of perfect graphene with H saturated 
edges. As an input for C-H and single/triple C-C 
bonds in the graphitic part and OPE, we use stan- 
dard literature values. For the interfacial contact, 
we take a single C-C bond with 1.39A. 

The choice of ribbons (width as shown in Fignj 
of single-layered graphite is guided by recent ex- 
perimental reports In those, self-assembled- 
monolayers of conjugated molecules were de- 
posited on a graphitic substrate. The molecular 
phenyl rings form a covalent bond to the carbon 
electrode in a similar fashion to that of Fig^ 
which ensures a low energy barrier at the in- 
terface. Our junctions is a first effort to model 
nanographite materials and may be considered 
as an idealisation with the following approxima- 
tions: a) the graphite material is well-ordered, b) 
the weak inter-plane interactions in graphite are 
insignificant, and c) there is weak inter-molecular 
interaction. These assumptions must be studied 
in detail separately, which is beyond our present 



scope. Later we discuss the validity of our results 
when other metals are considered as electrodes. 
Our motivation also stems from abiding interest 
in carbon-based nanojunctions (25126) . 

5. Results - Discussion 

We have calculated the conductance spectra for 
OPE wires consisting of up to five monomer units 
in a large energy window spanning the valence 
and conduction bands of poly-phenyl-ethynylene 
(PPE), using the formulation of SecEl In FigEl 
we plot for selected energies the logarithm of the 
transmission function for OPE molecular junc- 
tions of increasing length, represented by the 
number of monomer units N. The inset shows the 
local density of states at the carbon atom of the 
graphitic electrode which bonds to the oligomer. 
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Figure 2. Logarithm of the transmission function 
T{E) as function of the OPE length Na, for var- 
ious energies E. All energies are measured from 
Ep and lie in the band-gap of PPE. The inset 
shows the local DOS at the contact carbon atom 
of the graphitic electrode to the oligomer. 

There are two types of curves evident in Fig|2] 
One type, e.g., at E = — 0.75eV and 1.16eV, 
clearly indicates an exponential decrease of T{E) 
with length, from practically N = 1. This 
behaviour is robust for almost all energies in- 
side the gap formed by the renormalised HOMO- 
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LUMO energies of the embedded OPE. This gap 
is shghtly larger than the band-gap of PPE, al- 
lowing for the tunnelling behaviour. The other 
type of curve is approximated by Eq|21 after a 
minimum length; for E — 0.06eV at w 2 and 
for E = -0.09eV at TV « 3. The latter data 
show an anomalous length dependence with an 
initial increase of the conductance for A'' < 3. 
Such behaviour relates to unconventional MIGS 
|26| . These provide a way of molecular doping in 
short junctions that opens an additional resonant- 
tunnelling conduction channel. 

Unconventional MIGS are manifested when 
there is a large local DOS on the metal side. In- 
deed, the local DOS (inset of Fig[2Jl at the elec- 
trode interfacial contact of our molecular junction 
exhibits (van Hove) singularities. A most pro- 
nounced peak lies within the gap of PPE at the 
Fermi level of the system. This corresponds to a 
state located predominantly at the zigzag edges 
of the graphene ribbon and is responsible for 
the main departures of T{E) from Eq|21 Never- 
theless, the tunnel effect is still observed for quite 
small chain lengths and energies not in the im- 
mediate vicinity of local DOS peaks lOmeV), 
enabling the extraction of /3. This is a striking 
result since it extends the applicability of Eq|2] 
to ultra-short oligomer bridges and to a wide en- 
ergy spectrum not far-off to that responsible for 
resonant-tunnelling. 

We now turn our attention to the electronic 
structure of perfect PPE. Its complete band 
structure is calculated as described in SecOlvia 
EqdandlHI and is shown in FigOl The right panel 
corresponds to the conventional band structure. 
The spectrum of the imaginary part of complex 
k solutions, which come in conjugate pairs (±k), 
is plotted in the left panel. The complex-band 
structure shows the typical excursions which con- 
nect real bands between local extrema. 

An accurate complex-band structure calcula- 
tion is necessary for a comprehensive quantitative 
account of the MIGS [2H|. We have made ad- 
ditional calculations for alkane chains, graphene, 
and poly-para-phenylene (PPP). In all cases, our 
calculated band structures both in the real axis 
and in the complex plane, are in good agreement 
with those reported for either minimal or plane 



wave basis sets within DFT in the LDA |lll29l3()j . 

Moreover, as described in SecIHl MIGS are associ- 
ated with the complex k vectors for real energies 
and are exponentially decaying. Hence, we antic- 
ipate that states with the smallest |k| are domi- 
nant. Namely, of most importance for quantifying 
the damping factor in EqElis the, so-called, real- 
line connecting the top of the valence and bottom 
of the conduction bands, which occur at the bor- 
der of the Brillouin zone k — n/a. This line joins 
the energies Ey — — 1.12eV and Ec = 1.25eV, and 
is reproduced quite accurately within DFT-TB. 




Figure 3. Complex- and real- band structure of 
PPE in left and right panel, respectively. 

In Fig^ we isolate the positive branch of the 
Ec to Ey real-line and plot it on a larger scale. 
Together we plot the half of the damping fac- 
tor /3, which is numerically obtained by a least- 
squares fit to determine the slope of the lines in 
FigEl For energies where T{E) exhibits a non- 
monotonic length dependence, we have used only 
data for N > 2. The agreement is evidently re- 
markable for the entire energy spectrum, even at 
very small values of the tunnel decay parameter. 

There are hardly deviations from the relation 
f3{E) = 2 \k{E) \ a. K denotes that value of real- 
lines at E with the largest imaginary part in the 
renormalised HOMO-LUMO gap of the oligomer. 
Small irregularities appear only at those energies 
that the graphitic ribbons have a large local DOS. 
This marks the specific system and may not ap- 
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pear if strong disorder in graphene washes out 
the singularities. Since our results point to the 
intrinsic character of the damping factor for al- 
most any (3{E)N, we expect an even better agree- 
ment in the usual case of Al or Au electrodes. On 
the one hand, such broad-band metals usually ex- 
hibit seamless local DOS. On the other, larger 
interfacial barriers are formed in comparison to 
that of the C-C bond, due to the groups anchor- 
ing the metal. Small departures similar to those 
caused by the unconventional MIGS may occur 
if the anchor groups have localised states within 
the spanned energy window |11I29| . 
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Figure 4. The Ec to Ey real-line from FiglSl(solid) 
and the half of the damping factor as extracted 
from the slope of transmission function lines (see 
Fig|21) for all energies (dashed). 

Our results generally apply to oligomer bridges 
between metals in the off-resonant regime. They 
explain qualitatively the microscopic transport 
mechanism in experiments of ultra-short oligomer 
molecular junctions (4l5l(i| . and results obtained 
with non-equilibrium Green functions [31) . For 
tunnelling, an analysis ^B] similar to the Sim- 
mons model but including the spatial dis- 
creteness at these length-scales, may be used to 
derive current-voltage characteristics. Ep is es- 
timated by comparison of the extracted damping 
factor to the complex wavevector jT). 



Conversely, Ep is roughly located at the maxi- 
mum of the Ec to Ey real-line for metal-insulator 
(metal-semiconductor) junctions JT], which ap- 
proximately holds for the nanojunctions we stud- 
ied. Using this argument, we have calculated 
the damping factors for other polymers from the 
complex wavevector at the middle of their band- 
gap. For PPP we found P{E = 0) = 1.8, which 
is in quantitative agreement with the reported 
experimental I5I6I and theoretical 31 values in 
MMMs of oligo-para-phcnyl chains (OPP) up to 
four units between Au electrodes. However, one 
should take into account the possibility of non- 
planar intramolecular phcnyl-ring alignment in 
OPP. We shall report on such effects in forthcom- 
ing studies '33' and compare with other calcula- 
tions of the complex-band structure for twisted 
phenyl-ring geometries j29| . 

6. Concluding remarks 

We demonstrated that no tunnel-barrier col- 
lapse occurs even for ultra-short oligomer molecu- 
lar junctions, and for virtually all energies within 
the gap of the corresponding polymer, Ec — Ey. 
Contrary to expectations, there exists a unique 
damping factor /?(£'), even though /3{E)N 1 is 
not always valid. This is traced back to the larger 
value, in comparison to the band-gap Ec — Ey, 
of the renormahsed HOMO-LUMO gap of the 
oligomers in the junction. 

Moreover, we indicated the importance of cal- 
culating, and the need to understand and find 
ways to engineer, the complex-band structure 
of molecular wires. Our results establish that 
the damping factor /3{E) can be almost exactly 
mapped to the complex wavevector for all en- 
ergies within Ec ~ Ey. This substantiates that 
the complex-band structure provides a prolific 
method to determine the most crucial property 
of off-resonant electron transport across oligomer 
molecular junctions without performing a trans- 
port calculation. Improved independent first- 
principle methods may be used to increase the 
accuracy in evaluating (3. 

Our calculations on phenyl-ethynyl chains show 
that a single complex wavevector is adequate to 
explain the tunnelling regime but it cannot be ex- 
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eluded that a superposition of contributions from 
few close-by complex-A: may be needed to under- 
stand electron transport across other oligomers. 
Although there are strong indications that our re- 
sults apply in many cases, further investigations 
are clearly necessary to identify the typical be- 
haviour. In favour of the former are independent 
calculations on complex-band structures of vari- 
ous molecular wires, the rapid convergence of tun- 
nelling decay constants of short oligomers to the 
long wire limit (see also 131 35 ), and simple qual- 
itative arguments on the sensitivity of transport 
observables to small deviations in \k\ due to tun- 
nelling. 

The derived connection to the metal-induced 
gap states fortifies the analysis of the conductance 
properties of molecular junctions by using con- 
cepts developed in the theory of metal-insulator 
(metal-semiconductor) junctions. On the other 
hand, the tunnel properties of oligomers are of 
special relevance to molecular electronic devices 
and new circuit-designs. Finally, we note the 
conceptual similarity between electron transport 
in metal-molecule-metal junctions and the fun- 
damental chemical process of molecular electron 
transfer. Indeed, a formal analogy in a simple 
model relates the non-adiabatic electron transfer 
rate k^A in a donor-bridge-acceptor system to the 
Landauer conductance 0. The expression EqI21 
is reminiscent of the superexchange mechanism, 
and P coincides with McConnell's result |34I35| . 
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